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Abstract 

Delamination is a predominant failure mode in continuous fiber rein- 
forced laminated composite structures. One type of delamination is the 
transverse crack tip delamination which originates at the tip of transverse 
matrix cracks. An analytical model based on the sublaminate approach and 
fracture mechanics is developed in this paper to study the growth of such 
delaminations. Plane strain conditions are assumed and estimates are 
provided for the total strain energy release rate as well as the mode 1 and 
mode II contributions. The energy release rate estimates are used to 
predict critical delamination growth strains and stresses by assuming a 
critical energy release rate. These predictions are compared with experi- 
mental data on T30Q/934 Graphite Epoxy [±25/90 n ] s laminates in the range 
n=.5 to 8. A good agreement is demonstrated for the range of n where the 
experimental observations indicate transverse crack tip delamination to be 
the predominant failure mode. 

Introduction 

Fiber reinforced composites are now being used in a wide variety of 
engineering structures. The concept of directional strength and stiffness 
has been, for the most part, understood sufficiently to enable efficient 
load bearing designs. One of the current major issues in composite struc- 
tures is the understanding and prediction of damage modes and failure 
mechanisms. A thorough knowledge of the failure mechanisms is oound to 
lead to the design of efficient and durable structures. Failures in these 
materials often initiate in the form of matrix cracks or delaminations. 


Matrix cracks refer to intralaminar failures whereas delaminations refer to 


interlaminar failures. 

Matrix cracks usually occur within laminates where the fibers run at 
an angle to the primary load direction. Hence, such matrix cracks are also 
called transverse cracks. Based on the location and direction of growth, 
two distinct types of delamination can be discerned. These two types are 
called edge delamination and local or transverse crack tip delamination. 
Edge delaminations initiate at the load free edges of the structure whereas 
local delaminations start from a transverse matrix crack. In many cases, 
both types occur concurrently with varying levels of interaction. It has 
been observed in simple tension tests of uniform rectangular cross section 
specimen (Edge Delamination test) that delaminations initiate along the 
load free edges and propagate normal to the load direction. Transverse 
matrix cracks running parallel to the fibers have also been observed in off 
axis plies such as 90° plies. Such transverse cracks terminate where the 
ply orientation changes. Delaminations can originate at the interface 
where transverse cracks terminate. These delaminations, called transverse 
crack delaminations or local delaminations, grow normal to the transverse 
crack from which they originate. In the case of 90° plies, the growth 
direction is parallel to the load. 

The growth process of edge delaminations and local delaminations is 
often modelled using a fracture mechanics approach leading to the calcula- 
tion of a strain energy release rate. This is because the strain energy 
release rate can correlate delamination behavior from different loading 
conditions and can account for geometric dependencies. The strain energy 
release rate associated with a particular growth configuration is a measure 
of the driving force behind that failure mode. in combination with 
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appropriate failure criteria, the strain energy release rate provides a 

means of predicting the failure loads of the structure. 

Several methods are available in the literature for analyzing edge 

1-3 

del ami nations. These include finite element modelling , complex variable 
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stress potential approach , simple classical laminate theory based tech- 
nique and higher order laminate theory including shear deformations . 
Finite element models provide accurate solutions but involve intensive 
computational effort. Classical laminate theory (CLT) based techniques 
provide simple closed form solutions and are thus well suited for prelimi- 
nary design evaluation. Classical laminate theory based techniques provide 
only the total energy release rate, and thus in a mixed mode situation, 
there is insufficient information to completely assess the delamination 
growth tendency. A higher order laminate theory including shear deforma- 
tions has the ability to provide the individual contributions of the three 
fracture modes while retaining the simplicity of a closed form solution. A 
shear deformation model is available for edge delamination and has been 

/r 

shown to agree well with finite element predictions . 

Crossman and Wang^ have tested T300/934 Graphite epoxy [±25/90^3 s 
specimens in simple tension and reported a range of behavior including 

O 

transverse cracking, edge delamination and local delamination. O'Brien 
has presented classical laminate theory solutions for these specimen, 
demonstrating reasonabl e agreement in the case of edge del ami nation but 
with some discrepancies in the local delamination predictions. An empiri- 
cal finite element based combined edge and local delamination formulation 

g 

has also been proposed . Its predictions, however, do not fully explain 
the dependency of the critical strain on the number of 90° plies. 

In this paper, a shear deformation model is developed for the analysis 
of local delaminations originating from transverse cracks in 90° plies 


located in and around the specimen midplane. Plane strain conditions are 
assumed and thickness strain is neglected. Delaminations are assumed to 
grow from both ends of the transverse crack tip. The transverse crack is 
treated as a free boundary and the delamination is considered to be the 
crack whose growth behavior is to be modelled. The sublaminate ap- 
proach 1 ^’ 11 is used to model different regions of the specimen. The 

resulting boundary value problem is solved to obtain the interlaminar 
stresses, total strain energy release rate and energy release rate compo- 
nents. Critical local delamination growth loads are predicted for the 
[±25/90 n ] s specimen. 

Analytical Model 

The formulation is based on the sublaminate approach detailed in ref. 
10. A longitudinal section illustrating the geometry of a generic configu- 
ration is shown in fig. 1. The central region is assumed to be made of 90° 
plies with an isolated transverse crack in the middle. Oelaminations are 
assumed to grow from both ends of the transverse crack, and towards both 
ends as shown. From symmetry considerations, only one quarter of the 
configuration is modelled. The modelled portion is divided into four 
sublaminates as shown in fig. 2. The top surface (subiaminates 1 and 4) is 
stress free. In order to simplify the analysis, plane strain conditions 
are assumed and the thickness strain (e z ) is set to zero. The consequence 
of this combined with the fact that the w displacement is zero along the 
center line is that w is zero in sublaminates 1,2 and 3. Further, this 
approximation does not allow for the enforcement of boundary conditions on 
the shear stress resultants, leading to incorrect estimates of the inter- 
laminar normal stresses. The interlaminar shear stresses, however, are not 

fi in 

affected by this assumption ’ . The assumptions lead to considerable 

simplifications in the analysis. In spite of the simplifications, reliable 
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energy release rate components can be estimated based on the interlaminar 
shear stress distributions^’*^ 1 . 

A generic sublaminate is shown in fig. 3 along with the notations and 
sign conventions. The peel and interlaminar shear stresses are denoted by 
P and T respectively with t and b subscripts for the top and bottom surface 
respectively. The axial stress resultant, shear stress resultant and 
bending moment resultant are denoted by N, Q and M respectively. A summary 
of the governing equations is presented here for convenience. These are 
derived for a generic sublaminate using the principle of virtual work in 
Reference 12. 

The x and z displacements within the sublaminate are assumed to be of 
the form 

u(x,z)=U(x)+zp(z) (1) 

w(x,z)=W(x) . (2) 

Here U represents the axial midplane stretching and W is the transverse 
displacement. The shear deformation is recognized through the rotation 6. 
The origin of the coordinate axes for the sublaminates is taken at the 
delamination tip as shown in fig. 4. The equilibrium equations take the 
form 


N ,x + VV° 

(3) 

W p b=° 

(4) 

M x -Q+(h/2)(T t +T b )=0. 

(5) 


where h is the thickness of the sublaminate. The constitutive relations in 


terms of the force and moment resultants are 



where the A.., B. . and B-. are the classical laminate theory axial, cou- 
^ J ^ J ^ J 

pling and bending stiffnesses. The boundary variables to be prescribed at 
the sublaminate edges are 

N or U 
M or p 
Q or W. 

Additionally, at the interfaces between sublaminates, reciprocal traction 
and displacement matching boundary conditions have to specified. 

Solution Procedure 

A detailed solution is provided in the Appendix. A brief summary is 

provided here for convenience. The variables in sublaminates i and 2 are 

coupled by their reciprocal interlaminar stresses denoted T^ and P 1 and by 

displacement continuity at their common interface. Assuming exponential 

sx 

solutions for the axial force and bending moment resultants (N^=Ae , 
M 1 =Be sx etc.) leads to an eigen value problem involving the parameter s. 
The eigen values turn out to be 0 and two nonzero values (say s 1 and s 2 ) 
occurring in positive and negative pairs. Since the resultants maintain 
finite values as x tends to large negative values (left end of sublaminates 
1 and 2), the negative roots are dropped out of the solution. 

The following boundary conditions from the ends of the modelled region 
are enforced. 


n 2 (0)=0 

(9) 

Q 4 (a)=0 

(10) 

P 4 (a)=0 

(ID 

N l +N 2 = APPl i e d L°ad 

(12) 


Further, the following displacement matching conditions are applied. 

u 1 (x,-.5h 1 )=u 2 (x,.5h 2 ) 

u 1 (0)=u 4 (0) 
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(13) 

(14) 


u 2 ( 0 )=u 3 ( 0 ) 

(15) 

P 1 (0)=p 4 (0) 

(16) 


It should be noted that a ^ an< ( P 3 matching condition cannot be applied at 

this level of modeling since it would amount to specifying both W and 
6 12 

Q ’ . Consequently, there is a displacement discontinuity at the de- 
lamination tip. The effect of this will be discussed subsequently. To 

eliminate rigid body displacements, is set to zero at the left end. The 
following solutions can then be obtained for the resultants in sublaminates 
1 and 2 . 

N 1 =a 1 e s l x +a 2 e s 2 x +eA n(1) (17) 

N 2 =-a 1 e s l x -a 2 e s 2 x +eA n(2) (18) 

Ml= a i k i eSlX + a 2 k 2e S 2 x (19) 

M 2 =a 1 k 3 e s l x +a 2 k 4 e s 2 x ( 20 ) 

The interlaminar shear and peel stresses between sublaminates 1 and 2 can 
be obtained as 

T l= a l s ie SlX + a 2 s 2e S 2 X (21) 

P 1 =(k 1 +.5h 1 )(a 1 s 1 2 e s l x )+(k 2 +.5h 1 )(a 2 s 2 2 e s 2 x ) ( 22 ) 

In the above solutions, the k parameters are dependent on the eigen values 
and the stiffness of sublaminates 1 and 2 , the a parameters depend on the k 
parameters and the initial crack length a, and e is defined as 

c-o ( h ^+h 2 ) / ( A^^ ^ j+A^^ 2 ^ ) (23) 

where o is the applied uniform axial stress. Complete expressions for the 

eigen values and the a and k parameters can be found in the Appendix. 

Proceeding on to sublaminates 3 and 4, the following solutions can be 
written. 



V° 

(24) 


M 3 =<t>^sinh ^x+^^cosh w 3 x 

(25) 

where 

< * l 2'" a l* < '3 +a 2* < 4 ’ 

(26) 
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and 


^^“** 1^2 coth uj^ a 


0.5 


w 3 =(A 55(2) /D 11(2) ) 
N 4 =e(A ll(l) +A ll(2) ) 
M 4 =a l k l +a 2 k 2 


(27) 

(28) 

(29) 

(30) 


The corresponding displacement solutions are provided in the Appendix. 
The compliance of the specimen can be evaluated as 

C=2U 4 (a)/P (31) 

where P/2 is the load applied to the modelled section. The total energy 
release rate for the modelled section i.e. the total energy release rate 
per crack is then given by 

G t =P 2 /2 w (dC/da) (32) 

where w is the specimen width. Use of the previously described solutions 


leads to the following expression. 




(33) 


2w J (*4u(i) -Aii(i) + ^ii(2) + I' 2 

where the quantities and contain exponential terms dependent on the 
initial delamination length. Using the virtual crack closure technique, 
from the relative displacements in the cracked portion and the interlaminar 
stresses ahead of the crack tip, the mode I and mode II energy release rate 
contributions can be obtained. The mode III energy release rate is zero 
from the assumption of plane strain. The mode II energy release rate is 
given by 

n i: L [ t< ~ 

(34) 


G// = Um^ J T\(x — 6)t±u(x)dx 


where 6 is the virtual crack step size. The result of the limiting process 
is zero if there is no singularity in the stress field 10 . So, the limit is 
usually taken as the crack step size 6 tends to a small value, say a, based 
on the decay length or the length required to capture the essential fea- 
tures of the stress and displacement fields near the crack tip. The decay 
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In this study, the 


length is dependent on the eigen values s^ and s^. 
value of A has been set to 

Aq=. 25( 1/Sj + l/s 2 ) (35) 
since it reasonably fulfills the criterion given above. In a similar 
fashion, the mode I energy release rate can be obtained based on the normal 
stress (P) and the w displacements near the crack front. The normal (peel) 
stress estimate is inaccurate due to the absence of thickness strain. 
Hence, an alternate approach was used to estimate Gj, the mode I energy 
release rate. The total energy release rate for this problem is made up 
entirely of Gj and Gjj (Gjjj= 0). From an estimate of Gy and Gjj, an 
estimate for Gj can be obtained simply as 

G r G T ‘G II (36) 
The critical load for a given specimen, can then be evaluated based on an 
appropriate fracture law. This is illustrated in the following section. 

Results and Discussion 

The solutions derived in the previous section have been used to model 
the behavior of [±25/90 n ] s T300/934 Graphite Epoxy specimen for n values of 
.5, 1,2, 3, 4, 6, and 8. These correspond to the specimen tested by Crossman 
and Wang^. The specimen width and length were fixed at .0381 m and .015m 
respectively, as in the tests. The solutions were generated using a simple 
computer program based on the closed form expressions for the interlaminar 
stress and energy release rates. The applied load was set to 100 MPa, of 
the same order as in the tests. 

An example of the total energy release rate variation with the cracK 
length is presented in fig. 5. The asymptotic value of Gy is denoted by 
Gy 0 in the figure. It can be observed that after a certain crack length, 
the Gy is independent of the crack length. On the basis of curves like the 
one shown in fig. 5, the crack length was fixed at 10 ply thicknesses for 
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the remainder of the study. The dependence of the mode II contribution of 
the energy release rate on initial crack length (a) is depicted in fig. 6. 
Typical interlaminar shear and normal stress profiles are presented in 
figs. 7 and 8 respectively. The corresponding energy release rates have 
also been calculated and are presented in Table 1 and fig. 9. 

In order to evaluate the critical loads, an appropriate mixed mode 
fracture law has to be applied, based on the calculated energy release 
components. Since the calculated mode split shows only a small variation 
with n, the simple Griffith criterion G-j-=G-]- c has been used to scale the 

stresses to obtain the critical delamination growth stress (o c ) and strain 

2 

(e c ) values. The critical energy release rate Gy c was chosen as 415 J/m 
to obtain the critical stresses and strains listed in Table 1. This value 
of G^ c is larger than Gj c to account for the presence of mode II and the 
fact that G tt is about four times G T for the material system under 
consideration. The critical strains are plotted against n, the number of 
90° plies in fig. 10. The experimental results of ref. 7 and the predic- 
tions of refs. 8 and 9 are also presented in the figure for comparison. 
The predictions of the model developed in this paper are represented by the 
solid line while the experimental results are shown as filled squares. The 
classical laminate theory and finite element critical strain predictions of 
refs. 8 and 9 are represented by triangles with a connecting line and a 
dotted line respectively. 

In the experiments, the local delamination phenomenon was observed as 
the predominant failure mode only for the n=4,6 and 8 specimens. The shear 
deformation model presented in this paper provides good agreement with the 
experimental data in this range. For n<4, edge delamination either in the 
mid plane or in the 25/90 interface was observed in the tests. Hence, the 
predictions of the local delamination models in this region are not of 
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consequence as long as they do not predict critical loads lower than those 
predicted by edge delamination models. Thus, it can be seen that the shear 
deformation model predicts the observed behavior with reasonable accuracy 
and can be used in conjunction with an appropriate edge delamination model 
to predict critical loads accurately for the complete range of n values. 
The edge de lamination model presented in References 6 and 12 can be used 
for this purpose. However, a separate model is required to account for the 
mid-plane (Mode I) edge delamination behavior. 

Conclusions 

A shear deformation model has been developed to analyze local delami- 
nations growing from transverse cracks in 90° plies located around the mid 
plane of symmetric laminates. The predictions of the model agree reason- 
ably with experimental data from [±25/90 n J s T300/934 Graphite Epoxy lami- 
nates. The predicted behavior is such that, in combination with an edge 
delamination model, the critical loads can be predicted accurately in the 
range of n from .5 to 8. 
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Appendix A 

Sublaminate Analysis for Local Delaminations 

Interlaminar Stresses and Energy Release Rates 

A generic sublaminate is shown in figure 3 along with the notations and sign 
conventions. The interlaminar normal (peel) and shear stresses are denoted by P 
and T respectively with the t and b subscripts for the top and bottom surfaces 
respectively. The axial force resultant, shear force resultant and bending moment 
resultant are denoted by N, Q and M respectively. Plane strain conditions are 
assumed to prevail in the x — z plane and the thickness strain e zz is neglected. These 
assumptions lead to considerable simplification in the analysis. The displacements 
in the x and z directions are assumed to be of the form 

u = U(x) + z/3(x) (A.l) 

w = W(x) (A. 2) 

Here U represents the axial stretching and W is the transverse (thickness direction) 
displacement. This formulation recognizes shear deformation through the rotation 
(3. The equilibrium equations take the form 

= 0 


N, x + T t -T b 

Q,x + Pt ~ Pb 


o 


(A. 3) 
(A. 4) 


2 


M yX -Q + |(T, + T b ) = 0 (A. 5) 

where h is the thickness of the sublaminate. The constitutive equations in terms of 
the force and moment resultants are 


N — A n U tX + Bnfl tX 
Q = A 55 (ft + W <x ) 

M = BnU^ 4- D u ft yX 


(A.6) 

(A.7) 

(A.S) 


where A, B and D are the classical laminate theory axial, coupling and bending 
stiffnesses defined in the customary manner as 

An = J\C n dz 

h 

Bn = Cnzdz 


Dn = [\C n z'dz 
J 2 
r| 

A55 = J h C55 dz 

Here, the C s are the material moduli. For the case of plane strain in the x — z 
plane, the C s are defined as follows. 


(A.9) 


' 




f \ 

&XX 


C\\ Ci 3 0 


€-xx 

°zz 


Cl 3 C22 0 

< 

£zz > 

Txz 


1 

0 

0 

p 

Cn 


Txz 


The boundary quantities to be prescribed at the sublaminate edges are 

N or U 


M or ft 
Q or W 
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Farther, at the interfaces between sublaminates, appropriate reciprocal traction and 
displacement matching boundary conditions have to be used. 

The four sublaminates along with the loads acting on each axe shown in figure 4. 
Setting Pi and Ti as shown automatically satisfies the traction matching boundary 
condition at the 1-2 interface. From symmetry, we get w = 0 and zero shear stress 
along the bottom faces of sublaminates 2 and 3. This leads to w = 0 in sublaminates 
1,2 and 3. Thus, W has been prescribed in these sublaminates and the vertical shear 
force resultant Q cannot be prescribed at both ends of the sublaminates. Conse- 
quently, the calculated peel stress distribution will not be correct. In addition, at 
the 2-3 interface, the 0 s cannot be matched, since in these sublaminates, specifying 
0 is equivalent to specifying Q (through eq. A.7). Inspite of these simplifications, 
reliable energy release rate components can be estimated based on the interlaminar 
shear stress distributions. The mode I contribution can then be evaluated using the 
total energy release rate, which is not affected significantly by these simplifications. 

For the (±25/90 n ) s laminates under consideration, Bn is zero in all the four 
sublaminates. For sublaminates 1 and 2, the equilibrium equations and constitutive 
relationships can be written as 


£ 

H 

1 

II 

0 

(A.10) 

N 2tX + Ti = 

0 

(A. 11) 

Qi,x — P\ = 

0 

(A. 12) 

+ Pi ~ P 2 ~ 

0 

(A. 13) 

Mi, x + ^Ti-Qi = 

0 

(A. 14) 

M 2 , x + ^Ti-Q 2 = 

0 

(A. 15) 

H 

II 

f— * 


(A. 16) 
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N 2 — Au(2)C^2wc 

(A. 17) 

Q 1 = -^SS(l)/^l 

(A. IS) 

Qi = A 55 ( 2 ) 02 

(A.19) 

Mi = £„(,)&.* 

(A. 20) 

M 2 = D 11 (2)/?2 r r 

(A-21) 


The subscripts in brackets refer to the sublaminates to which the stiffness coefficients 
correspond. Equations A. 14, A. 15 and A. 12 can be rewritten in a modified form as 


M 1 , x + ^-N Ux = A 55(1) ft (A.22) 

M 2 ,X — ^N2,x = A 55 (2)/?2 (A. 23) 

Pi — Ql,x 

= M hxx + ^T ltX (A. 24) 

Matching the u displacement along the 1-2 interface implies 

Ul (“¥’*) = U2 (¥’ x ) 

or U t - = U 2 +!%/3 2 (A.25) 


Combining the equations to eliminate the displacement and interlaminar stress 
terms leads to the following homogeneous coupled system of ordinary differential 
equations. 





N 1 }X + N 2 ,x 

= 0 

(A. 26) 


Mi yXX + 

^Al,XX 

£>11(1) 

= 0 

(A. 27) 



t a ^- 

■^11(2) 

= 0 

(A.2S) 

iV, 

/z, M, 

n 2 

ho Mo 

= 0 

(A. 29) 

A n( i) 

2 -Dii(i) 

-4-11(2) 

2 £^11(2) 
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The solution is assumed of the form 


\ 

Ny 


' 4i ' 

n 2 

> = i 


Mi 


A 3 

m 2 

\ J 




>e' 


(A.30) 


Substitution of this solution into the governing equations results in the following 
system of algebraic equations. 

0 
0 




s 

0 


0 


0 

1 


L U(1) 




L n(2) 


2 _ ^ 55 ( 1 ) 
P H(1) 


0 

hi 1 


s 


2 D 


11 ( 1 ) 


2 _ 

•Dll(2) 

h"> 1 

TZ5 


11 ( 2 ) 



r 4i ' 


* \ 
0 

< 

4 2 

SX 

= 

0 


A 3 


0 




0 


(A.31) 


The corresponding eigenvalue problem has to be solved in order to obtain non trivial 
solutions. The eigenvalues turn out to be the roots of the following characteristic 
equation. 

s j-Bi-s 4 + B 2 s^ + B^ = 0 


(A.32) 


where 


By 

B 


1 

1 - 11 ( 2 ) 


1-11(1) 


1 

Di\(2) 


(t?) 

1 


2 = 


b 3 = 


1 


1 ..ffiglq. 

Ah( 2) -Dll(2) -^11(2) -^ll(l) 


4s5(l) , ^55(1) 

Unfit -Diini 77 


1 


^ll(l) 2- / ll(2) 
+ 


. - - ^5£ il + 1 - f > 55 l?J + j;55(2) 1 

-4-11(1) •C'll(l) 4-11(1) 1^11(2) B U (2) D n(i) 


§ 


1 4-55(1) 4-55(2) _ 1 4 55 (1) A 55 (2) 

4h( 2) -On(i) -Dll(2) 4-11(1) -Dll(l) £*11(2) 


For the material system and ply stacking sequence considered, B 2 > 4BiB 3 . Hence, 
the roots can be written as 


s = 0, ± 


— B 2 ± \J B 2 — 4 ByB 3 
2TT X 


(A. 33) 
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Only the zero and positive roots of eq. A. 33 are considered as they give exponentially 
decaying solutions, leading to finite values for the resultants at the sublaminate ends. 
Hence, the solution for N\ can be written as 


Ni = Oje* 1 * + a 2 e ,7X + oq 


Using this in eq. A.26 yields 


N 2 = —a ie 4ir — a 2 e'** + a 2 


(A. 34) 


(A.35) 


Substituting N x and N 2 in eqs. A. 27 and A. 28 provides the solutions for the bending 
moments as 


Mi = aikie SlX + a 2 k 2 e S7X 
hd 2 — aik 3 e six -)- o, 2 k^Q S3X 


The k parameters in the above solutions are defined as follows. 

¥4 


ki = 


'll(l) 


277 ^ _5 1 


k 2 = —Z- 


, 2 


fcs = 


TmT) $2 

¥4 

SS£21 _ 2 


^11(2) 
k < = -X 


¥4 


k'l 1 ( 2 ) 

If P is the applied force and w represents the specimen width, 


(A.36) 
(A. 37) 


(A.3S) 


(A. 39) 


(A. 40) 
(A. 41) 




(A. 42) 
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Using this in conjunction with eq. A.29 allows determination of the constants 
and a 2 . The following solutions for the stresses and the resultants can then be 
obtained. 


Nr 

_ .ill t - i P ^*ll(l) 

_ a t t + <.* +25 illl “+A 11(J) 

(A. 43) 

n 2 

- a r s »* n r a * x 1 P An(2) 

2u> A n (i) + A U ( 2 ) 

(A. 44) 

T\ 

= N hx 



= ajSje 4 * 1 + a 2 s 2 e SjI 

(A. 45) 

P i 

= M hxx + t$-T ltX 



= (fc + kt)ar$\e*' x + (k 2 + kj-)a 2 s\e* x 

(A. 46) 


The constitutive equations are used to write down the displacement solutions. The 
rigid body displacements of sublaminates 1 and 2 axe matched (in order to satisfy 
the displacement continuity condition) to obtain 


u ' ~ + & a 11(1) + (A - 47) 


u, = 


_ — r-^2 — e SJX + £■ 


L n(i)- S i 


•All(l) 5 2 


A 


11(1) 


+ ” 3 (A ' 4S) 


fix = + a 2 k 2 s 2 e a * x + + a 2 s 2 e s * x )} (A.49) 

fi 2 = a ^ — [ui^sSie* 11 + a 2 k4S 2 e 3iX + •^•(aiSie*' 1 -4* c^^e 521 )] (A. 50) 

•<*55(2) w 

The constants a x , a 2 and a 3 occurring in the solutions are determined using the 
boundary condtions. For sublaminate 3 the governing equations are 



= 0 

(A. 51) 

Qz,x + P 3 

= 0 

(A. 52) 

^3,x ~ Q3 

= 0 

(A. 53) 


£ 

cs 

II 

(A. 54) 

Ql — •d.55(2)^3 

(A. 55) 

= £>11 (2) /?3,x 

(A. 56) 


Matching U at the 2-3 interface and applying N 3 (a) = 0 gives 


N 3 = 0 


U 2 = U 2 (0) 


a i 


_£ 2 _ 


+ cl 3 


(A. 57) 
(A. 58) 
(A. 59) 


•51-^11(2) •S2-Ah(2) 

In order to solve for the bending moment, eqs. A. 53, A. 55 and A. 56 are combined 
to obtain 


M 3 ,xr - = 0 

■^ 11 ( 2 ) 

The solution of eq. A. 60 can be written as 

M 3 = <f> x sinhu; 3 :r -(- <j> 2 coshu> 3 x 
where the quantity u ; 3 is defined by 


(A. 60) 


(A.61) 


2 _ d. 55 (2) 
11 ( 2 ) 


^3 = XT 


(A. 62) 


(A. 63) 


Since the 0 matching conditon cannot be used at the 2-3 interface, the (remaining) 
boundary conditions are 

M 3 (a) = 0 
M 3 ( 0) = M 2 (0) 

The <f > s can be solved using the boundary conditions A. 63 as 

</f>2 = o, \k 3 + a 2 k^ 

h = 


— (j^cothu^a 


(A. 64) 
(A. 65) 
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The solution for sublaminate 3 can be completed by writing the following expres- 


sions. 


Qz = cosh o> 3 x -f <^ 2^3 sinhu^x 

(A. 66) 

fiz = [<^ 1^3 cosh u> 3 x + <^ 2^3 sinh u; 3 x] 

(A. 67) 

P 3 = tj— smhu; 3 x + <t> 2 coshu; 3 x] 
^1 1(2) 

(A. 68) 

The equilibrium equations for sublaminate 4 are 


£ 

H 

II 

o 

(A.69) 

Qa,x — 0 

(A. 70) 

M4 tX — Q\ = 0 

(A.71) 

The constitutive relations take the form 


N, — -4ll(l)C^4,r 

(A. 72) 

Qa — -4-55(l)(^4 + W 4tX ) 

(A. 73) 

Ma = -^ll(l)/?4,r 

(A. 74) 

Using eq. A.69 with the boundary condition N A (a) = yields 


II 

(A. 75) 

Similarly, using eq. A. 70 with Qi(a) = 0 results in 


<0 

Ox 

II 

0 

(A. 76) 

Matching M\ and M 4 at the 1-4 interface and using eq. A. 71 gives 


A^4 — d\}c j Cl 2^2 

(A. 77) 
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The U 4 displacement is obtained by integrating eq. A. 72 and using the displacement 
matching boundary condition U 4 ( 0 ) = £/i(0). 

u 4 = + + 33 ^) + “ 3 (A7S) 

Similarly, integrating eq. A. 74 and setting /? 4 (a) = 0 gives 


Pi = 7 ^ — [ai^i + a 2 k 2 )(x - a) (A.79) 

Using the solutions for Q 4 and /? 4 and the boundary condition W 4 (0) = 0 in eq. 
A. 73 yields the following solution for W A . 

W 4 = + a 2 k 2 }(^- - ax) (A. SO) 

In order to determine a l5 a 2 and 03, the following boundary conditions are used. 

JVi(0) = in 
MO) = MO) 

U\(—l + 0 .) = 0 


It is convenient to define the following parameters. 

*' = *!&<** + ¥> (*■“) 

*» - 1 (A - S2) 

* = 71^^+lt') ( A - S3 ) 

= rkt < A - S4 > 

&d = $3 — $1 ■+• ( 64 — 0 2 )a (A.S5) 

The nominal (far field) strain is given by 

„ _ P 1 


^ ^ 

Anjij 4- A U (2) 


(A.S6) 
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The a parameters are obtained as 

a i = ^ 11 ( 2 ) € \^ a(L (A. 87) 

a 2 = —^4.11(2) f (A.8S) 

u d 

«■ - < A - S9 > 

The specimen compliance C is defined as the ratio of specimen extension to applied 
load. This is obtained as 


r - 2U <( a ) 

= + “4 ( A£ >°) 

The total energy release rate associated with the crack (delamination) growth under 
a constant load P is given by 


Gf = 


_ P 2 dC 


tw da 


(A. 91) 


Using the compliance expression from eq. A. 90 in eq. A. 91 yields the following 
expression for Gj- 


r - P 2 

Gt ~s- 


(•An(i) in(i) + A U (2) 


(A. 92) 


where 


T — 1 -^- 11 ( 2 ) ^ 2^3 — f 1 — e _ 1 — e / A go) 

+ An(2) Aii(i) V s ! S 2 / 


T 1 ^n(2)(^3 + ^q)e- J>(< - a) --^ 1 +^Q)e- ,2(/ ~ a) 

2 ■4n( 1 ) + Au(2) 


(A. 94) 


The individual fracture mode contributions to the energy release rate can be cal- 
culated using the virtual crack closure method, based on the interlaminar stresses 
and displacements in the vicinity of the crack tip. From the assumed plane strain 
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condition, the mode III contribution is zero (G/j/ = 0). The mode II energy re- 
lease rate, G//, is calculated using the virtual crack closure technique while G/ is 
evaluated using 

G/ = G'/ — G// (A. 95) 

G// is calculated from the interlaminar shear stress and relative sliding displacement 
as 

G// = Hm J !Z\(z — 5)Au(i) dx (A. 96) 

In the absence of a singularity in the stress field, the result of the limiting process 
leads to the trivial result G// = 0. Hence, the limit is calculated as 6 tends to 
some finite value, say A. The value of A is chosen depending on the decay length 
associated with the problem i.e. the length within which the presence of the crack 
significantly alters the specimen response in comparison with the corresponding far 
field values. Evidently, the decay length in this problem is dependent on the eigen- 
values Sj and s?- The following value of A has been chosen in order to reasonably 
fulfil the decay length criterion. 



(A. 97) 


The relative sliding displacement An is based only on the difference U 4 — U 3 so that 
the kinematic condition of zero relative displacement at the crack tip is fulfilled. 
This also simplifies the calculations. The mode II energy release rate component is 
obtained as 

Gu = (A.9S) 

where J 3 and I 4 are defined as 


h = 


An(i) 


+ a ~ ) (ft ■ + §?) M 1 - - «-‘)] 


(A. 99) 
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Si A — 1 + e ,lA , s 9 A — 1 + e S2A \ + -^- 11 ( 2 ) 

a ' + ^ V in(„ 


(A. 100) 


Transverse Crack Spacing 


Shear Deformation Model 

The model presented so far has dealt with delaminations growing from a trans- 
verse crack. The same model can be modified to predict the spacing of these trans- 
verse cracks. In order to accomplish this, the delamination effect has to be isolated 
from the model. This can be achieved approximately by letting the crack length 
a tend to zero. This yields an approximation since the boundary conditions axe 
not accounted for properly by this limiting process. To get an accurate shear de- 
formation model, we consider only sublaminates 1 and 2 and apply the following 
boundary conditions for sublaminate 2. 

N 2 ( 0) = 0 (A. 101) 

M 2 (0) = 0 (A. 102) 

Using these boundary conditions in eqs. A. 37 and A. 44 yields two equations in a x 
and a 2 which can be solved to obtain 


_ k d P An( 2 ) 

1 fc 4 — k 3 2w An(i) + An(2) 

(A. 103) 

k 3 P Au(2) 

' x2 — k 4 2w An(i) + An( 2 ) 

(A. 104) 


The interlaminar shear stress can now be obtained using eq. A. 45. The saturation 
crack spacing corresponds to the distance from the crack where the broken plies 
regain their uniform stress/strain state i.e. where the interlaminar shear stress has 
decayed down to its far field (uniform) value. Practically, this distance is calculated 
by looking for the x where the interlaminar shear stress is some small fraction (say 
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.001) of its maximum value. The maximum shear stress evidently occurs at x = 0 


and is given by 


T i(mar) , 

1 — Gl$l + CL 2$7 


(A. 105) 


The crack spacing A can then be determined by solving the following transcendental 
equation. 


aT.Sie JlA + a?s?e ,2A 
GjSl + C2$2 


0.001 


(A. 106) 


Membrane Model 

A simpler model can be used to estimate the saturation spacing of the transverse 
cracks. This model treats the sublaminates as membranes i.e. the bending effects 
are ignored. The equilibrium equations for a generic membrane sublaminate are 


N s + Tt-T b = 0 (A. 107) 

%(T t -T b )-Q = 0 (A.10S) 

i 

The constitutive equations take the form 

N = ( A ' 109) 

Q = A ss /? (A. 110) 

The displacements are assumed to be of the following form. 

u = U(x) + z/3(x) (A. Ill) 

w = 0 (A. 112) 

The following governing equations can now be written 

Ni yX — Ti = 0 (A. 113) 

N 2 , i + T 1 = 0 (A. 114) 
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1 fT l -Q 1 = 0 

(A. 115) 

^Tx-g 2 = 0 

(A. 116) 

Nx = jiU\' X 

(A.117) 

Ni = 72^2,1 

(A.11S) 

Q 1 = ■^■55(1)01 

(A. 119) 

Q 2 = A 5 s(2)02 

(A. 120) 


(A. 121) 

where the 7s are defined as 


»2 

71 = Al, W uSS) 

(A. 122) 

72 = All(2 > iS 

(A. 123) 

Eqs. A. 113 and A. 115 can be combined as 


Qi = 

(A. 124) 

Using eqs. A. 119 and A. 117 in this leads to 


0i — _ 2 l ~a — ~l\Ui, xx 

Z -^55(1) 

(A. 125) 

Following a similar procedure for 02 yields 


& = hrr' u '** 

Z -^55(2) 

(A. 126) 

Using these two relations in eq. A. 121 leads to 


Mfc) Ar, u '~ = u '-&) tEsP'- 

(A. 127) 

Combining eqs. A. 113, A. 114, A.117 and A. 118 gives 


llUi iXX + 72 U 2 , XX = 0 

(A.12S) 
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Substituting this into eq. A. 127 results in 


U\' XX 






*■ 55 ( 1 ) \ * J ^- 55 ( 2 ) 

The characteristic roots of this differential equation are 


\yXXXT 4" i^Ul^xx 0 


s = 0, 0, ± 


n 


7i +72 


7i72 


3^ + (%) Z 1 


55 ( 2 ) 


(A. 129) 


(A. 130) 


The solution for U\ can then be written as 


U i — Aie six + A 2 X + Az (A. 131) 

where the As are arbitrary constants to be determined from the boundary condi- 
tions. The root S\ is the positive root such that a decaying solution is obtained in 
the negative x region. For the special case of 5 n (i) = Bn ( 2 ) = 0, the nonzero roots 
can be written in a simpler form as 

1 

— |_ -2 .. 

^ 55 ( 1 ) ^ 55 ( 2 ) 

The interlaminar shear stress can be obtained as follows. 


(A. 132) 


s 2 = jjylgl) 


Ti = N hx 

= 71 ^ 1 , XX 

= 7iAiSje iia: (A. 133) 

The maximum shear stress is 

Ti max) = ll A 1 sl (A. 134) 

Then, the saturation crack spacing A corresponds to 


e S| A = 0.001 


(A. 135) 
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Shear Lag Model 

This model allows for a nonlinear displacement field through the thickness of the 
sublaminate. Its fundamental assumption is that the shear deformation neglected 
in the classical theory of bending can be estimated using the shear stress. The 
sublaminate axial force equilibrium condition can be written as 

N x + (T t -T b ) = 0 (A. 136) 


The axial stress is assumed to be uniform and is given by 



The shear stress is estimated as follows 


(A. 137) 


&xz y z — 

' -A r .x 

= ~H~ 
= 

rt 


(A.13S) 


This can be integrated to obtain 

a xz = lL^ z + Ii^3± (A. 139) 

Neglecting transverse displacement, the axial displacement can be obtained by in- 
tegrating the shear strain, which in turn is obtained from the shear stress. 


u , = 


tT 

0-55 

1 


(Tt - T 6 )f + 21+^ 


u 


= U(x) + ^\(T t -T b )£ + (T t + T b )z 


(A. 140) 
(A. 141) 


where U(x ) is the mid-plane axial displacement. This displacement expression can 
be used to obtain an improved axial stress estimate as follows. 


® XT 


G 1 1 ^ ,X 


IS 


= c, 


n 


U' x 2CsS^ 1 ~ ^' x ~h "** ^ xZ 


(A. 142) 


The corresponding axial stress resultant can be written as 



(A. 143) 


The governing equations for the sublaminate are thus eqs. A. 136 (equilibrium), 
A*141 (displacement field) and A. 143 (constitutive relationship). Using these to 
model sublanunates 1 and 2 results in the following governing equations. 


Ni tX — T\ = 0 
-^2,x + T 2 = 0 

Ni = (7x1(1) 

N2 = Cn( 2 ) 


hiUi ' x - 

h? 


\h 2 U 2iX + ^cl (2 Ti, x 
Ui = Ul + 2 + 

U 2 = U2 + 2ck^[ Tl £ + TlZ . 

Displacement continuity at the 1-2 interface implies 

= u 2 (x,-^) 
or U 2 = U t - ^ 

Equation A. 146 can be rewritten as 

TT — Ni y hi m 

U '* - + 54C^ t -.' 

Combining eqs. A. 147, A. 151 and A. 152 results in 


AT (~* f h 2 N} h 2 Ty. x 

A2 " cn(2) l 


V 


hi j_ -J 12 . 


55(1) ^55(2) J 


(A. 144) 
(A. 145) 
(A. 146) 

(A. 147) 
(A.14S) 
(A. 149) 


(A. 150) 
(A. 151) 


(A. 152) 


(A. 153) 
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But from eqs. A. 144 and A. 145, we have 


N 2tX = -Ti = —N\ iX 


(A. 154) 


Using this in the differentiated form of eq. A. 153 leads to 


TuC 


2'~'11(2) 


+ 


ThC 


n(i)J 


Ni 


_ 1 


+ 


1 ^ 55 ( 1 ) ^ 55 ( 2 ) 


N- 


l, XXX 


(A. 155) 


The nonzero characteristic roots of this equation are given by 


2 _ o ( C ^ + ^2^11(2) A 

VAlCii(i)y V ^2^11(2) / V^2^55(l) + ^1^55(2 )/ 


(A. 156) 


This is the same as in the membrane model except for the factor 3 which is 4 in the 
membrane model. This difference is related to the fact that the axial displacement 
distribution through the thickness is parabolic in the shear lag model and linear in 
the membrane model. The crack spacing A for the shear lag model is determined 
as in the case of the membrane model but using the modified characteristic root. 



Table 1 Summary of Results 


number of 
90° plies 

Gt 

J/m 2 

Gii/Gj 

o- c 

MPa 

e c 

% 

1/2 . 

2.404 

0.276 

1313.9 

1.6747 

1 

6.752 

0.275 

784.0 

1.1685 

2 

22.849 

0.267 

426.2 

0.8058 

3 

51.049 

0.261 

| 

285.1 

0.6427 

4 

93.603 

0.256 

210.6 

0.5444 

6 

228.871 

0.250 

134.7 

0.4264 

8 

440.065 

0.247 

97.1 

0.3555 
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Fig. 2 Modelled Region and Sublaminate Scheme 
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